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Abstract. 

Whereas Lie had linearized scalar second order ordinary differential equations (ODEs) 
by point transformations and later Chern had extended to the third order by using 
contact transformation, till recently no work had been done for higher order or sys¬ 
tems of ODEs. Lie had found a unique class defined by the number of infinitesimal 
symmetry generators but the more general ODEs were not so classified. Recently 
classifications of higher order and systems of ODEs were provided. In this paper 
we relate contact symmetries of scalar ODEs with point symmetries of reduced sys¬ 
tems. We define new type of transformations that build up this relation and obtain 
equivalence classes of scalar third order ODEs linearizable via these transformations. 
Four equivalence classes of such equations are seen to exist. 


1 Introduction 


Symmetries of ODEs have acquired considerable interest over the last few decades. 
The earlier work had dealt with second order linear ODEs. The number of point 
symmetries a scalar second order ODE can have is exactly one of 0,1, 2, 3, or 8 [IS] , 
Lie mmmm discovered that if a second order scalar ODE possesses 8 point 
symmetries, then it can be converted to a linear equation ( linearizable ) via invertible 
point transformations. He proved that a scalar second order linearizable ODE is at 
most cubic in the first derivative. The coefficients of linearizable ODEs satisfy four 
constraints. These constraints were reduced to two by Tresse [23]. 

Mahomed and Leach ca showed that all linearizable scalar ODEs of order n with 
(n > 3) can be put into three equivalence classes with n + 1, n + 2 and n + 4 Lie point 
symmetries. Chern DEI used contact transformations to reduce scalar third order 
ODEs to y'" = 0 and y"' + y = 0. Grebot pj |6J addressed the linearization problem 
of scalar third order ODEs via restricted point transformations. Neut and Petitot 
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[ 20 ], and independently Ibragimov and Meleshko (IM) [7J [8] used the original Lie 
procedure of point transformations to linearize scalar third order ODEs to the form 
y"'[x ) + ko(x)y(x) = 0. They determined the linearizability criteria as well as the 
procedure for the construction of linearizing transformation for scalar third order 
ODEs. IM also used contact transformations to linearize scalar third order ODEs. 
The linearization problem for scalar fourth order ODEs gets more complicated and 
was dealt with by Ibragimov, Meleshko and Suksern (IMS) [ 3 ( 22 ]- They obtained 
the linearization criteria both by point and contact transformations. Meleshko [T9] 
also provided a simple algorithm for reducing third order autonomous ODEs to Lie 
linearizable second order ODEs and was extended to scalar fourth order autonomous 
ODEs [3]. The class of equations linearizable by Meleshko algorithm are not included 
in IM and IMS classes. 

As such there are three equivalence classes for scalar third order ODEs linearizable 
via point transformations but no classification for these equations linearizable via 
contact transformations was provided. Indeed IM [8] obtained the linearizable form 
via contact transformations for a scalar third order ODE but they neither discussed 
its contact symmetries nor provided the classification for it. The subject of the 
present paper is to investigate the contact symmetries of linearizable scalar ODEs. 
In fact we reduce the scalar third order ODE to a system of two second order 
ODEs and relate contact symmetries of the scalar ODE with point symmetries of 
the reduced system. In doing so we define a new class of transformations that are 
more general than contact transformations and perform group classification of scalar 
third order ODEs linearizable via these transformations. 

The plan of the paper is as follows. In the subsequent section we provide some 
preliminaries. In the third section we define generalized contact transformations. 
The canonical form of linearizble scalar third order ODEs is obtained and group 
classification is performed by using this form. The last section contains a summary 
discussion of the results. 


2 Preliminaries 


A differential equation (DE) is said to be invariant ( symmetrical ) under a point 
transformation 


t = ip(x,y), u — if(x, y), 

if it is form invariant under that transformation. Instead of using finite transforma¬ 
tions we can use infinitesimal symmetry generators 


X = £(z,2/) 


d_ 

dx 


+ v(x,y) 


d_ 

dy 
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To apply the above symmetry generator on n th order ODEs, we need to extend or 
prolong it. The prolonged generator np to order n is given by 


x (n) = £(®, y )^~ + v (x, y)^ + v {1) (®, y, y')-^ + v {2 \x, y , y ', y")^ 


dy 


dy' 

+ ... + y { - n \x,y,y',...,y { - n) ) 


dy" 

d 

dydd ’ 


where rj^ is the extension coefficient of and is given by 



dr f n 1} y (n) d i 

dx dx 


If we have a system of m ODEs of order n i.e., m dependent variables y = Vi, (i = 
1,2,... ,m), then 


X (,,) = £{x,y i )^ + 7y(x,y i )-^ + ri 1 \x,y i ,y-)-^+ y ( l 2 \x,yi,y , i ,yi) 1 ^jj 


dy , 


d y[ 


dy- 


■ (n) / / (n) \ 9 

+ --- + vl {x,yi,Vi,...,y\ )—7 




(n) 


(1) 


Now will be the symmetry generator or simply the symmetry of the system of 
m ODEs of order n 



p / in (n— 1 )\ 

fi{x,y i -,y i ,y i ,...,yl 


(i = 1,2,..., m), 


( 2 ) 


if it satisfies the relations 

X (n) [^ n) - fi(x, yi, y[, y",..., yi n ~ 1] )} = 0, (i = 1,2,. 
which imply that 

7 /$ n) = X (n - 1 )/*, (z = 

which are called symmetry conditions for the system of ODEs 


,m), 


( 3 ) 


Two equations are said to be equivalent if there exists a point transformation which 
maps one equation into the other. All first order scalar ODEs are equivalent. How¬ 
ever this does not hold for second or higher order ODEs. Nevertheless all linear 
scalar second order ODEs are equivalent [16] and can be transformed to y" = 0. For 
scalar n th order linear ODEs with n > 3 we have a theorem due to Laguerre mm 
which says that any linear homogeneous n th order scalar ODE 

n— 1 

y ( n) + ^ ki(x)y {l) = 0 , n > 3, (4) 

i =0 

can be transformed by a point transformation to an equation of the form 

71 — 3 

y<"> + y>(;i)y (i > = 0. (5) 

i =0 
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Eq (J5J) is called the Laguerre canonical form of the linear scalar ODEs flU) . Hence 
there are three equivalence classes arising from the Laguerre canonical form for the 
scalar third order ODEs. These three classes are: y"'(x ) = 0, y"'[x ) + k 0 y(x ) = 0 
and y"'(x) + ko(x)y(x) = 0, which have 7, 5 and 4 point symmetries respectively. 

We need to generalize contact transformations. Hence we shall first review the 
definition of these transformations. A transformation 

t = <p(x,y,z), u = ip(x,y,z ), v = w(x,y,z), (6) 

where z = y\ is called a contact transformation if it satisfies the contact condition, 


3 Generalized contact transformations 


Consider an n th order scalar ODE with (n > 3) 

</"> = f(x,y, (7) 

We substitute y' — z in the above ODE. This reduces the above scalar ODE to the 
following system of two ODEs of order n — 1 

y(n~ i) — 

z^ =f(x,y,z-z',z",...,z^). (8) 

A point transformation 

t = <p(x,y,z), u = if(x,y,z), v = u(x,y : z), (9) 

for the system (J8]) corresponds to a generalized contact transformation for the scalar 
n th order ODE (J7J) with z = y'. The transformation (J9]) is actually the contact 
transformation without the contact condition. 

The system of equations (JHJ) is linearizable if there is a point transformation (!9|) 
which reduces the system (JHJ) to a linear system of ODEs. 


3.1 Group classification 

Consider the general form of a linear scalar third order ODE 

y'" = S(x) + a (x)y + a(x)y' + /3(x)y". (10) 
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Defining y' — z we reduce the order to 2 and double the dimensions, i.e. 

y" = d, 

z" = S(x) + <r(x)y + a(x)z + fi(x)z' . (11) 

Here the variable z is actually the derivative of the variable y. For the purpose of the 
group classification we replace the variable z by y' and compare the corresponding 
system of equations (fTTl) with the second canonical form of the linear system of 
second order ODEs [24J given by 

y" = ki(x)y' + k 2 (x)z', 

z" = k 3 {x)y' + k±(x)z'. (12) 

Since (1T2T) only depends explicitly on y' and z' so it gives 5(x) =0 and a{x) = 0. 
The system (TTT1) is now of the form 

y" = 

z" = a(x)y' + /3(x)z'. (13) 

This is the reduced form of the systems of second order ODEs that is obtained from 
a scalar linear, third order ODE. We now perform the group classification on the 
above system. For this we suppose 

j., , <9 , . d . . d fp. . d 

> = £{x,y,z)— + rn{x,y, z)— + y 2 {x,y, z)— + y\ J (x,y,z)— 

+y 2 } (x,y,z)-^ +Vi\x,y,z)-^ +4 2) {x,y,z)-^j, (14) 

be the symmetry generator for (fT3l) . then symmetry conditions (J3J) read as 

i i? = x ( 1 | (z'), 

^* 2) = X a) (a(x)y + I3(x)z‘), (15) 

/ m /o\ 

where r/, and 772 J are second extension coefficients of the generator (fT4l) and are 


given by 

hi — yi,xx T y (2hi ,xy i\,xx ) 4- 2z T/\ t xz T y (hi ,y ^£,x 3y y 2z ^z) 

4 ~Z (hi ,2 y £,z) T y (jll,yy 2^xy 2 Z ^,,yz) T 2 y Z (hi, yz £,,xz) 

+^ /2 (hl ,22 - y't,,zz) - y' 3 {,,yy , (16) 

h2 — y2,xx + Z (2hl,a :z £, xx ) T 2 y T]2,xy T Z (j]2,z 2 £ x 2 y t^ y 3 Z £ i2 ) 

T y (h2 ,y Z £,j/) T % (h2 ,22 22 Z ^,yz) 4“ 2 y Z ( J]2,yz £,m/) 

+h' 2 (h2,22 - z'iyy) - Z' 3 t,zz- (17) 

The symmetry conditions (fl5l) give the following system of determining PDEs 

£,2/2/ = £,yz = £,zz = 0, y 1, 2Z - £,* = 0, T}l }XX - T) 2 ,x = 0, ( 18 ) 

hi,2/2/ — 2^x2/ — a £,z = 0, h2,2/2/ ~~ a £,2,/ = 0, ‘^’> 12 ,zz ~ 2£,xz ~ £,2/ — 2 / 3 £,2 = 0, ( 19 ) 

2hi,2/z - 2e,xz - Ky - Piz = 0, 2 rn,yz ~ K*y ~ Ply ~ 2 <z = 0, (20) 

2hl ,xy ~ r) 2 ,y - £,,xx + arji >z = 0, 2 t] 2 ,xy + Ctr) 2 ,z - Py 2 ,y - atr)i,y - £a, ^ - a^ x = 0, (21) 

2 hl ,xz 4 "" hi ,y h2,z £,a: 4 “ dhl,z = 0 , 2 ?^ 2 ,xz «hl,z 4 “ h2,2/ £,xx £P,x d£,x = 65 ( 22 ) 

m,xx - ayi,x - Pm,x = o. (23) 
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The system of PDEs (fT8|) gives the following solution 


£ = yai(x) + za 2 (x) + a 3 (x), 

( 24 ) 

m = 2 z2 °2 + za4,(x, y) + a 5 {x, y), 

( 25 ) 

\ 

% = 2 ^ ^2,# ^^4,# H - 0-5,x ^6 (2/? 

( 26 ) 


where a*, (i — 1, 2 ,..., 6) are arbitrary functions of their arguments. 


We now assume f3{x) to be zero, nonzero constant and an arbitrary function of x 
and consider the following cases. 


Case I / 3 (x) — 0 

The system of PDEs (1191) — fl 23 [) in this case takes the form 

Vl ,yy 2£ ; X y ,yy ®£,y = If 2t)2,zz 2£ xz y — 0, (27) 

2m,VZ - Kxz - Ky = 2m,yz - K*y ~ 2 <z = 0, (28) 

2lil,xy - m,y - £,xx + Otm.z = 0, 2m,xy + OiT]2,z ~ am,y ~ ~ = 0, (29) 

2m,xz “t m, y y?,z £,x = 0; 2m,xz ayi,z fr m,y £ ,xx = o, (30) 

V 2 ,xx - am,x = o. (31) 

The above system of PDEs is solved for different values of ot(x). 


Case 1.1 Both a and (5 are zero 

In this case the system of PDEs (l27ll — (l3Tf reduces to 


Vl,yy 2£ <X y — 0 ) 92 ,yy — 0 , 2t)2 ,zz 2£ xz £ ,y — 0 , 2 ? Jl,yz 2£ xz “1£ y — 0 , 

2r)2,yz 2^ X y — 0, 2?]l,xy V2,y £,xx — 0) 2r\2, X y — 0) 

2ri\, xz + T)\ y r)2,z £,x — 0) 2ri2,xz T r\ 2 ,y £,xx — If d2,xx — 0- 

Solving the above system yields the following 15 Lie point symmetries: 


X, = 

X4 = x 
X 7 = x 


_ 9 _ 

dx’ 

d_ 

W 

d 


X 2 = 


dy 


X, = 


dz 


X 5 = z 


1 2 3 

dx 2" dy ’ 


d 1 2 9 d 

fry’ 6_ 2 X ~frj +X !fr’ 

d 1 d 1 

X ’ =x d^ + r%' X8 = 2“ 

1 2 d d d 

Xl1 = w- + xy-^~ + y-z-, 

2 dx dy dz 

v 2 d 1, 2 \ d 9 

X12 = x — + -(x z + xy)— + xz 

dy 


d 


X 10 = i^xz + y)-^-, 

2 dy 


dx ' 2' 

. 

-xz- 


‘dz 


-1 d 1 .9 1 2 d 

x “ = T“aS + P“ + 2!/) 5S + 2 2 a? 

„ 1 .9 1 . 9 

X.4 = ( 5 m + ») s - i *(x 2 -2 ! ,)-, 

v /1 2 > 9 1 2 2 2\ d ,1 2 J 

- 9) 5; + ( 2 12 - yz) ai‘ 


d 


dy *dz’ 


(32) 

(33) 

(34) 

(35) 

(36) 

(37) 

(38) 

(39) 
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Case 1.2 a = a 0 ^ 0, (3 = 0 

The system of PDEs (l2?j) — (j3Tj) in this case yields a 15—dimensional Lie algebra. The 
first three operators are X l5 X 2 , X 3 , given by (132]) . while the remaining 12 operators 
are 


v 9 d 

Yi = y— +z—, 
oy oz 

Y 2 = Y.3 = e ^|- + 

oy oz oy oz 


d 


'Oz 
d 1, 


Yt = v ii + es k + ^ y2a ° +i:2) i' 

Y5 = *s + 5 (#2 “° + * 2) i + “'> 2! 4- 

o n 

Y 6 = e-^"— - 

oy oz 

o o a 

Y 7 = — + V^e^ x y— + a 0 e^ x y—, 

Y 8 = - ^T 0 e-^ x y^ + a 0 e~^ x y^ 

{zy/a^ - a 0 y)e^ x d mx d 

Yg = —--—-— + (zVdo - «oy)e v 0 

yao dy dz 

(2^o7 + ao/de - ^* d y ~^ x 9 

Yio = —--—- — - (Zy/ao - a 0 y)e V°“> x 

y/oo oy Oz 

Y (v/aio y — z)ed“° x d (a 0 y 2 — z 2 )e^° x d ( a 0 y 2 — z 2 )e^° x d 

11 = 7^0 a® + V^o d^ + 2 d? 

(y/aoy + z)e^d“° x d (a 0 y 2 — z 2 )e~^o x d (a 0 y 2 — z 2 )e~^° x d 


\/c*o 


da; 


2 v / oii dy 


+ 


dz' 


(40) 

(41) 

(42) 

(43) 

(44) 

(45) 

(46) 

(47) 

(48) 

(49) 

(50) 


Case 1.3.1 a = (x ± c) m , m ^ 2, or e x , /3 — 0 
This case produces a 5—dimensional Lie algebra. The first three operators are 
X 2 ,X 3 , Yi given by (l32j) and (l40lh 


Case 1.3.2 a = (x ± c) 2 , (3 = 0 

In this case we obtain a 6—dimensional Lie algebra. 


Case II j3{x) ^ 0 

The following subcases arise: 

Case II. 1 a = 0, j3 = /3 0 ^ 0 

The system of PDEs (ITDl) — fl23|) simplifies to 

Vl ,yy ‘l^xy = 0j ?? 2 , 2 /y = 0? 27/2 ,zz ‘l^xz = 0j 

27/1,150 2 £,,xz 2£,j5 Po^,z = 0) 2 7/2,1/2: do Ly = 0) 

27/1,312/ ^72,2/ = 0? 27/2,a;2/ do??2,2/ = V2,xx do^?2,oi = 0> 

2?7l, 12 T ^7l,2/ ^?2,2 £,x 4“ 3oVl,z = 0? 27 J2,xz 4“ ??2,2/ £,xx do£,a; — 0? 

which produces a 7—dimensional Lie algebra. The first four of these operators are 
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X 1 ,X 2 ,X 3 , Y 4 given by (1321) and (1401) while the remaining three are 



(-A) V + Z )j!-, y 4 

dy 


eP° x d 
Po dy 


, J3 0 x 

+ dz ’ 


Case II. 2.1 a = a 0 = /3 0 ^ 0 

In this subcase the system of PDEs ffTO l) — (1231) gives the following set of solution 

£ = ci, m = c 2 y + c 3 z + c 4 + c 5 e aix + c 6 e a2X , 
m = c 5 a \e aiX + c e a 2 e a2X + c 3 (y + z)a 0 + c 2 z + c 7 , 


where 

«i = ^(«o + \A*o + 4a 0 ), «2 = ^(a 0 - y/«o + 4 «o), 

and q, (z = 1, 2,..., 7) are arbitrary constants. This yields a 7—dimensional Lie 
algebra with X 1; X 2 , X 3 and Y x given by (1321) and (1401) . The extra three operators 
are 


v d d 

Y 2 = 2— + a 0 [y + z ) —, 
ay dz 


Y, = 


d d 

+ Y 4 

ay oz 


d_ 

dy 


e a2 X —+a 2 e a2X 


d_ 

dz' 


Case II. 2.2 a = a o ^ 0 and /? = /5o 7^ 0 with a 0 ^ fio 

This case produces the following set of solution for the PDEs (TlDli — fl23l) 

£ = ci, g i = c 2 y + c 3 z + c 4 + c b e PlX + c & e p2X , 

12 = c 5 Pie^ lX + c 6 / 3 2 e fex + c 3 a 0 ?/ + c 3 / d 0 ^ + c 2 z + c 7 , 

where 

Pi = ^(A> + + 4a 0 ), and /3 2 = ^(do - + 4a 0 ). 

From above we get a 7—dimensional Lie algebras. The first four operators are given 
by (J32D and (|40|). 


Case II.3 (3 = /3 0 , a(x ) = cad 71 , (x + c) m , m = 1,2 

This case produces a 5—dimensional Lie algebra with first three generators X 2 , X 3 , Y 4 . 


Case II.4-1 a(x) = (cx ± d ) m , f3(x) = (cx ± d ) m , m = 1, 2 
Here we find a algebra with 6 Lie point symmetries. 

Case II.4-2 a(x) = P{x) = a^ _1 ),a^ -2 ) 

This case lies in the above case. 


The system of PDEs (TT8T) — f|23l) provides us four equivalence classes with 5, 6, 7 
and 15 symmetries. These Lie point symmetries correspond to generalized contact 
symmetries for the scalar ODE (fTOjh Thus we have the following theorem. 

Theorem 3.1. If a scalar third order ODE is linearizable via generalized contact 
transformation (J9]) then it has one of 5, 6, 7 and 15 generators of generalized contact 
transformations. 
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4 Conclusions 


Though Leach and Mahomed ca had shown that there are three equivalence classes 
of third order scalar ODEs linearizable via point transformations, no work on sym¬ 
metry group classification of these equations linearizable via contact transforma¬ 
tions was done. In fact IM [8] got the necessary form of a scalar third order ODE 
linearizable via contact transformations but there was no attempt to address the 
classification problem with their methods. In this paper we have found a connec¬ 
tion between (generalized) contact transformations of system of order n with point 
transformations of system of order n — 1. By defining the first derivative of de¬ 
pendent variables to be new variables we reduce the order of a system from n to 
n — 1 and increase its dimension form m to 2 m. Point transformations for the lower 
order system correspond to generalized contact transformations for the higher order 
system. We obtained the canonical form of scalar third order ODEs linearizable via 
generalized contact transformations. This canonical form gave us four equivalence 
classes for scalar third order ODEs depending on the number of infinitesimal gen¬ 
erators. As such, a scalar third order linear ODE has three classes with 4, 5 and 
10 contact symmetries. The maximal symmetry class with 15 generalized contact 
symmetries corresponds to the maximal symmetry class of contact symmetries. Also 
of great interest is the study of correspondence between the other classes of general¬ 
ized contact symmetries and those of contact symmetries. Here we obtained group 
classification of a scalar third order ODE by reducing it to a system of two second 
order ODEs. If the reduced system of ODEs is linearizable then it can be solved by 
using geometric linearization ED- 

We can reduce a scalar fourth order ODE to a system of two third order ODEs and 
following the same procedure of group classification we can find generalized contact 
symmetries of the scalar ODE. Similarly we can take it to higher order ODEs to get 
equivalence classes of these equations by simply reducing scalar ODEs to systems of 
ODEs. One could use the given procedure to find the equivalence classes of systems 
of higher order ODEs. We can reduce a scalar fourth order ODE to a system of 
two second order ODEs by double reduction. Similarly any system of ODEs of 
order n > 3 can be reduced in steps to a system of second order ODEs to use the 
benefit of geometric linearization. This procedure can also be employed in steps 
to relate higher order transformations to point transformations [4]. There is much 
work that needs to be done in this direction. The equations linearizble via these 
transformations may form a new class of ODEs that do not fall into IM or IMS or 
Meleshko classes of linearizable ODEs. 
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